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Abstract 



An analysis of semileptonic decays of B mesons with the emission of a single 
soft pion is presented in the framework of the heavy-quark limit using an 
effective Lagrangian which implements chiral and heavy-quark symmetries. 
The analysis is performed at leading order of the chiral and inverse heavy 
mass expansions. In addition to the ground state heavy mesons some of 
their resonances are included. The estimates of the various effective coupling 
constants and form factors needed in the analysis are obtained using a chiral 
quark model. As the main result, a clear indication is found that the 0^ 
and 1^ resonances substantially affect the decay mode with a D* in the final 
state, and a less dramatic effect is also noticed in the D mode. An analysis 
of the decay spectrum in the D^*^ — vr squared invariant mass is carried out, 
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showing the main effects of including the resonances. The obtained rates show 
promising prospects for studies of soft pion emission in semileptonic 5-meson 
decays in a i?-meson factory where, modulo experimental cuts, about 10^ such 
decays in the D meson mode and 10^ in the D* mode could be observed per 
year. 



2 



Typeset using REVT^ 



I. INTRODUCTION 



Semileptonic B ^ D^*^ decays with emission of a single pion may soon be established at 
CLEO or ARGUS, as well as at the planned i?-meson factory. These decays, denoted B14 in 
the rest of this article, complete the list of this category of decays which includes and 
Di^. There are fundamental differences among these three decays, which make their study 
very interesting. In particular, Ki^ decays can be studied within chiral perturbation theory 
(xPT) over the whole final-state phase space as the resulting pions are soft Di4 decays 
101 are much harder to study because, with the exception of a small fraction of phase space, 
the final state involves light mesons with relatively large energies. This makes the use of an 
effective theory not viable. 

decays offer new theoretical possibilities. As a whole, they are difficult to predict 
because the kinematic domain of the daughter pion ranges from the soft limit, to be properly 
defined later, to a high energy limit. However, it is possible to restrict the study to the soft- 
pion limit, where one can make use of the powerful constraints resulting from heavy quark 
spin-fiavor symmetry and chiral symmetry. 

In this work we study the soft-pion domain in Bg^ decays, for which the effective chi- 
ral Lagrangian approach combined with the inverse heavy mass expansion p| provides a 
consistent theoretical framework. In a strict sense, it turns out that the soft-pion domain 
represents about 80% of the Bg^ rate in the D mode and about 10% in the D* mode. 
Here the large fraction in the D-mode is mostly due to the inclusion of the cascade decay 
B — > £uD* — > IvDtx. Since, according to a rough estimate, the total branching ratio for the 
5^4 decay is about 1% (0.2 %) in the D (D*) mode, it seems that experimental access to the 
soft-pion domain in the foreseable future, such as at a i? meson factory, is clearly possible. 

Perhaps the most compelling motivation for studying Bi^ decays is the overall current 
status of the semileptonic decays of B mesons. The measured inclusive semileptonic branch- 
ing ratio is 11.1±0.3 %, while the sum of measured exclusive semileptonic branching fractions 
is significantly less than this. In particular, the elastic modes B —>■ Dev and B D*eu 
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account for only about 60 % of the total semileptonic branching fraction |^ (6 ^ m modes 
are expected to be suppressed by |Vb„P/|VbcP ~ 0.01). Understanding the so-called in- 
elastic modes, such as those we discuss here, is therefore crucial to resolving the apparent 
discrepancies among the measurements. 

Another interesting aspect of -B£4-decays is that they may give an indication of resonance 
effects (especially D meson resonances). At present, the only well established resonances 
are the two P-wave objects Di and 1^2, with = 1"*", 2"*", respectively. States that may 
contribute significantly to B^^ decays, and hence which may be observed in such decays, 
include the remaining lowest-lying P-wave states (with = 0+, 1+), two of the lowest lying 
D-wave states (J^ = 1^, 2~), and the first radially excited S-wave states {JP = 0~, It 
turns out that the well established Di and D2 states do not play a role in the present 
analysis, while only the S-wave radially excited D mesons offer any opportunity for direct 
discovery using the Bi^ decays, as they are the only ones with a small enough width to show 
up as a resonance feature in the decay spectrum. 

From the theoretical standpoint, the soft-pion limit is particularly interesting, as chiral 
symmetry and spin-flavor symmetry determine to a large extent the different decay ampli- 
tudes, in terms of a few low-energy constants and universal form factors. These form factors 
are associated with the matrix elements of the charged 6 — c electroweak current between 
the relevant heavy meson states, and the low-energy constants determine the amplitudes of 
the strong interaction transitions mediated by the soft pion. 

Another area of interest in decays is their contribution to p, the slope of the Isgur- 
Wise function which, in turn, impacts on the extraction of \ Vcb\ from data. Bjorken, Dunietz 
and Taron have obtained a sum rule that relates the slope of the Isgur-Wise function 
for the elastic decays B —>■ Diu to the form factors that describe the inelastic semileptonic 
decays. The decays that we consider here provide the leading resonant and non-resonant 
contribution to this quantity. 

These decays have been analyzed in the framework we explore by a number of authors. 
B — > Duel' has been treated by Lee, Lu and Wise [0], and by Cheng et al. while Lee 
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0] and Cheng et al. have examined B D*Tiev. In these analyses, only the ground 
state mesons, the D, D*, B and B* were included, which amounts to keeping only those 
contributions wich are leading at the zero recoil point v ■ v' = 1, v and v' being the four 
velocities of the B and D mesons, respectively. In addition, we note that only Cheng et al. 
went on to estimate branching fractions and evaluate the differential decay spectra. 

In this work we include resonances which contribute at leading order in the chiral expan- 
sion, and neglect those which correspond to radially excited states (except the two resonances 
0~ and 1~ which we have to keep as shown later). Since resonance contributions vanish at 
zero recoil in the infinite heavy quark mass limit, their effects can only be observed away 
from that point, and as we will demonstrate, they alter substantially the decay rates. 

The predictions that arise in our analysis rely heavily on our ability to obtain good 
estimates of the aforementioned low-energy constants and universal form factors. In this 



work we use the chiral quark model [T^ convoluted with wave functions obtained in a 
simple model of heavy mesons to give such estimates. It is our hope that this procedure 
leads to reasonable results. 

The experimental status of -8^4 decays is hazy. ARGUS has analyzed the process B — > 
D**£^u as background to the semileptonic decays B Diu and B D*ii>. The resonances 
included in their analysis were the four P-wave states alluded to above, as well as the two 



radially-excited S-wave states. They report 63± 15±6 possible candidates [|TT|. This result 
is based on studying the invariant mass distribution of the D*tt combinations that result 
from the decay of the D**. The corresponding branching ratios are BR{B^ D**£^u) = 
(2.7 ± 0.5 ± 0.5)%, when their results are fitted to the model of Isgur, Scora, Grinstein and 
Wise (ISGW) [in, and BR{B^ D**tu) = (2.3 ± 0.6 ± 0.4)% when fitted to the model 
of Ball, Hussain, Korner and Thompson (BHKT) [|13[. This result implies that the resonant 
contribution to decays is significant. In the case of the decay B D-Kiu, the D* 
provides the dominant contribution, largely because of its proximity to the Dtt threshold. 
As far as we know, no other experimental group has published numbers for semileptonic 
decay rates of B mesons to excited D mesons. 



This article is organized as follows. In section II we review the effective theory resulting 
from spin-flavor and chiral symmetries. Section III presents the calculation of the effective 
coupling constants and form factors appearing in the effective theory, while we present the 
analysis of the decay amplitudes and differential widths in section IV. Section V is devoted 
to the results and discussions. A number of calculational details are relegated to three 
appendices. 



II. EFFECTIVE THEORY 

In this section we briefly review the effective theory which incorporates simultaneously 
spin-flavor and chiral symmetry and describes the interactions between soft pions and 
heavy mesons. Within the framework of the heavy quark effective theory (HQET) [|T^ , a 
hadron of total spin J consists of a light component (the brown muck) with spin j, and the 
spin-1/2 heavy quark, with J = j ± 1/2. For a given j, there are therefore two mesons, and 
these are degenerate members of a spin-flavor multiplet, at leading order in HQET. In the 
rest of this article, we denote the multiplets by the of the two states. For example, for 

= 1/2^, we have the multiplet (0^, 1^). 

For reasons that we will outline later, the only multiplets of interest in this work are 
(0~,1~), (0+,l+), (1",2-), and (0", 1")'. Here, (0", 1")' denotes the first radially excited 
version of the ground state (0^, 1~) multiplet. In order to formulate the effective theory, 
it is very convenient to introduce superfields associated with each multiplet [|1^]. These 



superfields provide a natural way of realizing the spin-flavor symmetry. At leading order 
in the inverse heavy mass expansion one associates one such superfield with each possible 
four velocity f^. This is because in the large-mass limit of the heavy quark, a velocity 



superselection rule sets in 
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The superfield assigned to the ground-state heavy meson multiplet (0,1) with velocity 

^- = ^ (-75P + 7'%*), (2.1) 
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where P and V* {v^V* = 0) are the fields associated with the pseudoscalar and vector 
partners, respectively. These fields contain annihilation operators only, and are obtained 
from the relativistic fields as 



P(x) = ^/M e-'^''-'' <l>W(x), 

V;{x) = e"^^^^-^ ^^^^\x), (2.2) 

where the label (+) refers to the positive frequency modes of the relativistic field, and M is 
the meson mass. 

The spin-symmetry transformation law is 

7i _ exp S^„) 7i _ , 



Si = ^e^''y,Jl]^-^^. (2.3) 

where e^, k = 1,2,3, are space-like vectors orthogonal to the four-velocity. 7i_ = 7o7il7o 
transforms contravariantly to 7i_. 

In a similar manner, it is straightforward to define the superfields associated with the 
excited states Jl^. In our case, the states of interest are the (0^^, 1+), (1^, 2^) and (0^, 1^)' 



multiplets, which are described by the superfields 

1+ ^ 



\ V ^ L O J 

1+ ^ 



^- = ^ (-75i^o-+7.i^r-)' (2-4) 

respectively. All the tensors are transverse to the four-velocity, traceless and symmetric. 
The transformations of these superfields under spin symmetry operations are implemented 
in exactly the same manner as in the case of 7i_. 

The chiral transformation law of the superfields is easily determined by following the 
well known Coleman- Wess-Zumino procedure to implement non-linear realizations of non- 
Abelian symmetries. All multiplets are isodoublets (we do not include the s-quark in our 



analysis), so that the transformation law under an arbitrary chiral rotation belonging to 
SU{2)l^SU{2)r is 

h{L,R,u)n, (2.5) 

where Ti. is any spin symmetry multiplet, and h{L,R,u) is an SU{2) matrix which results 
from solving the system of equations 

Lu = uh{L, R, u), 

Ru^ = u'^h{L,R,u). (2.6) 

Here, L (R) is an SUl{2) {SUr{2)) transformation and u is given in terms of the Goldstone 
modes (pions) as 

u{x) = exp (^-^n(x)^, 

U{x) = lf-T', Fo = 93MeV. (2.7) 

The transformation ( p.5|) is like a gauge transformation, the x-dependence entering via the 
Goldstone boson field. In order to build an effective Lagrangian which is chirally invariant, 
a covariant derivative is thus required, and is 

T, = Tl = ^ iud.u^ + u^d.u) = ^ [n, d,U] + .... (2.8) 

Another fundamental element in the construction of the effective Lagrangian is the pseu- 
dovector, 

i 1 

cj^ = - {ud^u^ - u^d^u) = j^dfj,U + . . . , (2.9) 

which transforms homogeneously under chiral transformations. 

Since the rest mass of the different heavy mesons is removed according to equations 
analogous to eqn. (|2.2|), acting on the respective superfields is proportional to the 
residual momentum carried by the superfield, which in the present analysis is of the order 
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of the pion momentum. This imphes that V^^ counts as a quantity of 0{p) in the chiral 
expansion. Obviously, a;^ is of the same order. 

Throughout this work we consider only the leading terms in the inverse heavy mass 
expansion and in the chiral expansion, and neglect the effects of chiral symmetry breaking 
due to the light quark masses (they enter only via the pion mass when the final-state phase 
space is considered). As we point out later, we have to include the effects of hyperfine 
splitting in the ground-state mesons, even though this splitting represents a departure from 
the spin-flavor symmetry. These modifications are only kinematic. 

To leading order in xPT the amplitudes for the Bu decays are proportional to a single 
power of the pion momentum. This places a restriction on the angular momentum quantum 
numbers of the excited states that can be considered in this analysis. The excited states 
that contribute to the B14 decay amplitudes at this order can be identified by examining 
their decays, via soft-pion-emission, to the ground state heavy mesons. Consider an excited 
state with spin J — j ± 1/2, where j is the spin of the light component of the meson (the 
brown muck). Let us define the integer k = j — 1/2. The soft-pion decay amplitude of such 
a state to the ground state supermultiplet is of O (^p'^^ if the parity of the excited meson is 
(— l)'^, or of O (p'^'^^^ if the parity is (—1)*^+-'^. In the case where k — and the parity of the 
resonance is positive, as is the case of the (0^, 1+) multiplet, the amplitude is proportional 
to p - V and is of order p as well. Using these 'rules', one finds that only the states with the 
quantum numbers mentioned above contribute at leading order in xPT. 

The effective theory is given by an effective Lagrangian which explicitly displays spin- 
flavor and chiral symmetry. The 0{p) strong interaction part is 



= Tr (d^Ud^U'') - \b Fl Tr {M{U + C/^)) + O(p'), 
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£x" = -i^^^TiD V ^-.j + . . . 

C^- = -^^/.TrD (:^'_ V'' + . . . . (2.10) 

C/(x) = li^ix) and Ttd is the trace over Dirac indices. The interaction terms involving o;^ 
have not been displayed in the Lagrangians of the resonant states, as they are not needed 
in this work. The only interaction terms in we need are those which give transitions 
between the resonances and the ground state mesons via a single soft pion emission. These 
are characterized by three low energy constants ai, a2 and a^-, and are 

>C™ = CKi Ttd 

+ aa TrD {fi!. uj^H-^^n^ . (2.11) 

The vertices resulting from needed in this work are displayed in Appendix A. 

The range of validity of the soft-pion limit is estimated from the invariant product of 
the pion momentum and the four velocity. As long as this product is smaller than some 
scale A^, which is of the order of 0.5 GeV, the application of the soft-pion limit should 
be appropriate. In this limit, this product has to remain smaller than the mass splittings 
between the neglected resonances and the ground state mesons, thus giving an estimate of 
the value of A^. Since two velocities appear, namely the velocities of the parent B meson 
and the daughter D meson, the soft-pion limit requires that the invariant products of the 
pion momentum with both velocities must be smaller than A^^ 

Another set of essential ingredients are the matrix elements of the electroweak charged 
current C7^(l — 75)6. Since this current is an isosinglet, it does not have direct couplings to 
pions at leading chiral order, and its matrix elements are easy to parametrize in the effective 
theory. Denoting the superfields corresponding to D and B mesons and their resonances 
respectively by X> and B, the matrix elements of the charged current are obtained from the 
effective current operators 

(a) J^((0-, 1-) ^ (0-, 1-)) = i{v) TVd (^-(^')7m(1 - , 
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(b) J^((0+, 1+) ^ (0-, 1-)) = p,{v) [Tro {V^W)^,{1 - 75)i3-(^)) 

+ Tro {V-{v')^,{l-^,)B+{v)) , 

(c) J^((l-, 2-) ^ (0-, 1-)) = p2{y) [tvTrc {v'+{v')l,{^ " 75)i3-(t;)) 

+ t;;TrD (P_K)7m(1-75)S^M)], 
(d) J,((0-, 1-)' ^ (0-, 1-)) = Tro (:n'_K)7,(l - lr>)B.{v) 

+ P_K)7M(l-75)S-(t;)). (2.12) 



Here, v {v') is the four- velocity of the B [D) meson, and v = v-v' . ^{u) is the Isgur-Wise form 
factor, which is normahzed to be unity at zero recoil (z/ = 1) if one ignores QCD corrections 
and higher orders in the inverse heavy mass expansion. The other form factors, namely 
Pi{y) and P2{i^)i which determine the transition between the resonances and ground 
state mesons via the charged current, are not constrained by symmetry at zero recoil. The 
currents however do vanish at zero recoil due to kinematic factors. In fact, if a resonance is 
characterized by a given value of the previously defined parameter fc, the current of interest 
is suppressed by a kinematic factor of the form {y — f%j {v — v')^^...{v — v')^,,. Note that 
no suppression of this form appears for the current of eqn. ( p.l2| (d)). However, heavy quark 
symmetry and orthogonality together imply that ^'■^''(z^) has to vanish at zero recoil, at least 
as {u — 1). The expressions for the currents of eqn. ( p.l2| ) are given in Appendix B. 

We note that by choosing to work to leading order in p, we have also placed an implicit 
restriction on the powers of u — 1 that appear. Since the largest value of k to be considered 
is unity, we find that the amplitudes for the decays are proportional to at most a single 
power of f — v', and the differential decay width will contain terms with at most two powers 
of z/- 1. 

We conclude this section by making a final comment on the states we include in our 
analysis. As outlined above, working to order p has severely restricted the states we can 
include, at least as far as their angular momentum quantum numbers go. However, there is 
no restriction on their 'radial' quantum numbers. Our self-imposed restriction of excluding 
any radially excited states (with the exception of the radially excited (0~, 1~)' multiplet) 
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is motivated by two related factors. These radially excited states are expected to be quite 
a bit more massive than their non-radially-excited counterparts. Thus, we expect little 
contribution from such states, provided we do not venture too far from the non- recoil point. 
Furthermore, the propagators of such states are expected to lead to a further suppression of 
any contribution, as in the strict soft pion limit these states will be far off their mass shell. 



III. LOW ENERGY CONSTANTS AND FORM FACTORS 

The soft-pion limit of the -8^4 decays is determined in terms of four low energy constants: 
g, ai, (X2 and 03; four universal form factors: ^(z^), ^^^^{y)-, P\{y) and P2{y)\ and mass 
differences between the resonances and the ground state mesons: bm\ = Mq+ — Mq- , 6m2 = 
M2- — Mq-, and Sm^ = Mg/- — Mq-, and the total decay widths of the resonances. In writing 
this form, we are treating the states in each excited multiplet as degenerate with each other. 
However, in dealing with the contribution from the ground-state doublet (0^,1~), it is 
imperative that we include the mass difference SniQ = Mi- — Mq-, as this plays a profound 
role on the outcome of our analysis. We begin by formulating a simple model of the heavy 
mesons, and using the wave functions obtained from this model to calculate the quantities 
we need. 

To estimate the masses and obtain wave functions, we use a version of the Godfrey-Isgur 
model UTS). In this version, one of the quark masses is set to infinity. In addition, we do 



not expand the wave function of a state in a harmonic-oscillator basis, but instead choose 
it to be that of a single harmonic-oscillator state with the appropriate quantum numbers. 
The oscillator parameter of each wave function is obtained in one of two ways. In the first 
method, we perform a variational calculation, and the values obtained in this way are listed 
in table |I[ Also listed in this table are the mass differences between the excited states and 
the ground states. We discuss the second method of obtaining the value of the gaussian 
parameter below. 

The low energy constants appropriate for soft-pion emission are estimated in a chiral 
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quark model. In this model, pions couple only to the light constituent quarks of a heavy 
hadron, via the Lagrangian 

H = iq-i^Wf,q - rhgqq + gA{^)qit,l->uJ^q, 
Vf,q = df,q + iTf,q, (3.1) 

where and cu^ were defined previously, rhg is the constituent quark mass, and the con- 
stituent quark field q transforms under chiral rotations as g ^ h{L, R, u)q. The axial 
coupling of the quark, (7^(0), is assumed to be unity. Arguments favoring this value for the 



axial coupling constant have been given in [|T^]. A non-relativistic expansion of the interac 



tion term of this Lagrangian is performed, and the resulting non-relativistic interaction term 
is convoluted with wave functions obtained from the model described previously. The low 
energy constants obtained in this way are also listed in table |I[ More details of this model 
will be presented elsewhere. 

One of the results of this analysis is that some of the low energy constants vanish in the 
limit when the pion energy in the vertex is taken to zero. For the cases where this chiral 
suppression occurs, we define the low energy constants as corresponding to the energy of 
the pion in the decay of an on-shell heavy resonance to an on-shell heavy ground-state. We 
expect that this procedure will furnish only rough estimates of these couplings. 

We also use the chiral quark model to estimate the total and partial widths of the excited 
states relevant to our analysis. Phase space limits all decays to pions or r^'s, both of which 
can be described in terms of chiral dynamics. It turns out that decays to ?7's play only a 
small role, and then only for the (1~,2~) multiplet (F(i- 2-)^(o-,i-)77 ~ 4 MeV). 

We can compare the results we obtain for the partial and total widths of these states 



with the calculation of Godfrey and Kokoski ||2^ , for instance. Unfortunately, we have only 



a single pair of states in common with that calculation, namely the (0''", 1'*') multiplet. The 
large total widths we obtain for these states are consistent with the widths obtained in pO 



The low energy constants ai, a2 and are related to the partial widths for the 
resonance decays into ground state mesons via single pion emission by 
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2 "tt' 



ro+ — ri+ 

f 1- = f 2- 

To-' = Ti-/ 



3a? 






Mn + 6mi 


a\ 


Md 




Md + 5m2 




Md 







pI 



(3.2) 



It is interesting to note that the value g = 0.5 obtained here is a direct result of the 
assumption that g\{0) = 1 and is independent of the wave function used. 

The total widths of these states are similar to the partial widths obtained in this fashion, 
with one exception. The total width of the 2~) resonances is dominated by their decay 
into the (1+, 2+) resonances, with a resulting total width of 405 MeV. 

The form factors ^, pi and p2 are also obtained using these wave functions. They 
are extracted from the overlap of the wave function of the ground state with the boosted 
wave function of the appropriate excited state. The boost we use is a Galilean boost, which 
means that we are neglecting relativistic effects, as well as effects that arise from Wigner 
rotations. The explicit forms we obtain for these form factors are 

A2 



exp 
1 A 



- 1) 



4/?2 

4/32 

' 2/5/5' 



1 

5/2 



v/2/3 V/9' + /5'V 
1 /AV ( 2(3(3' 



exp 

exp 

\ 7/2 



A2_ 

4/52 



A^ 



2 (/?2 + (3' 



'/2\ 



exp 



A2 



2 (/32 + /?'2) 



(3^3) 



In these expressions, (5 and (5' are the harmonic oscillator parameters of the ground and 
excited states, respectively. A is defined by writing the mass of the ground state as M(o-,i-) = 
mg + A. In the second of 

eqn. (^.31), we have set (3 = (3' to ensure orthogonality of the wave functions of the 
(0", 1") and (0-, 1")' multiplets. 
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The parameters obtained by the methods outhned above will be refered to as set I. We 
obtain a second set of values for (3 and jS' (and for all of the quantities we need, except the 
masses of the states) by first setting all the /5's to the same value, and then choosing this 
value so that it reproduces the experimentally measured slope of the Isgur-Wise function. 
The values of the parameters we obtain in this way are also shown in table |, and we will 
refer to this set of parameters as set II. 
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TABLES 



TABLE I. Quark-model parameters and low-energy constants used in this work, sets I and II. 



Multiplet 


(3 (GeV) 


M-M(o-,i-) (GeV) 


r (MeV) 


coupling constant 


(o-,i-) 


0.57 








0.50 


(0-,1-y 


0.57 


0.56 


191 


0.69 


(0+,l+) 


0.56 


0.39 


1040 


-1.43 


(l-,2-) 


0.51 


0.71 


405 


-0.14 



Multiplet 


13 (GeV) 


M-M(o-,i-) (GeV) 


r (MeV) 


coupling constant 


(0^,1-) 


0.5 








0.50 


(0-,1-y 


0.5 


0.56 


174 


0.66 


(0+,l+) 


0.5 


0.39 


756 


-1.22 


(l-,2-) 


0.5 


0.71 


408 


-0.215 
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IV. Bu DECAY AMPLITUDES AND DIFFERENTIAL WIDTHS 



A. B DttIu decay amplitude 



The decays we consider are D'^iun'^, B^ D^iun'^, B^ D^£u7i^, and 

B^ D^iun^, wliose amplitude magnitudes are in the ratios \/2 : 1 : -\/2 : 1. In what 
follows we give the results for the 7r° in the final state. The amplitude for this process has 
the general form 



(4.1) 



Here is the V-A charged leptonic current. For all practical purposes the lepton mass can be 
neglected (we do not consider decays into the r family) and the leptonic current is considered 
to be conserved. The hadronic part of the amplitude Q'^ receives non-resonant (NR) and 
resonant (R) contributions, illustrated in figures |] (a) and (b), respectively. Using the results 
of Appendices A and B, the evaluation of the Feynman diagrams is straightforward. The 
non-resonant portion is 



i v"v'^ e^a/sp + g^,p{l + i^) - Vf, v' 



+ 



i v'^v'^ e^ai3p + g^,pil + i^) -Vp v' 



(4.2) 



2{p-v' — 6m£)) + ie 

here 6m£) = mu* — rrio and 6m b = itlb* ~ itlb are the (positive) hyperfine splittings in the 
ground state multiplets, and 



(4.3) 



The resonant portion of fi^ is 



-TO 



p-v 



+ 



p-v 



(y.2 



2(—p-v — 6mi) 2{p-v' — 6mi] 



2{—p-v — 61712) 



+ g^^iu^ - I) - vl v^i2 + u)~3v' 
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+ 



2{p-v' — 6m2) 



Fo [—2p-v — orris 



+ ief,uapv''v"^{u - 1) + g^,{u^ - 1) 



i v^v'f^ e^app + g^p{l + v) - Vf, v' 



+ 



i v°'v''^ epapp + gpp{l + iy) -Vp v' 



(4.4) 



2p-v' — 6m3 

Here we denote 6m j = 5m j — iTj/2, where Tj is the total width of the resonance. 

As explained earlier, contributions from other resonances than the ones considered are 
suppressed either by higher powers of the soft-pion momentum, as is the case with the 
(l"*", 2"*") multiplet, by higher powers of (z/ — 1) in the small recoil domain, or by the fact 
that they are much heavier than the ground state mesons. 



B. B D*7r£u decay amplitude 



The amplitude for this process has the general form 



(4.5) 



where eD*u is the polarization vector of the D*. 

The non-resonant contributions to Q^'^ are obtained from the diagrams shown in figure 
|10| (a), which give 
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P" 2Fq^^ ' r "^p-v' + 5mr) + ie 



+ 



-{p-v + dms) + ie 
p-v' + ie 



9ua{V + V)p- Qp^Vy - QpyV^ + «e^„<^^(W + V 

i(^ppapv"^v'^ + gpp{l + v) - Vpv'^ 
The resonant contributions obtained from the diagrams in figure [10| (b) are given by 

- 9pu{^ - + v'pVy - iepuapv°'v''^ 



(4.6) 
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- 9ua{v + v')f,{l' - 1) + ^V^v'^v'^ 



- '^9nuvU^ - 1) + gt^aVuii" - 1) + ^e^,aau{v - 't;')"(l + 1^) 



+ 3p^ 



2{p-v' — Srh2) + ie 



- Va [v'^{2 + iy)- 



V6gpi,{v - 1) - vsVpv'^ + itf.sapVpv'^v''^ 
ie^^aapV^v'^^iiy - 1) + g^a{i^'^ - 1) 



Pv 



+ p 
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p_%M 



2{—p-v — Sfhs 



^2{p-v' — 61713) 



(4.7) 



2{p-v' — 5777.3) 

Here, 0^^ results from the numerator of the spin-2 propagator in the heavy mass limit, and 
is 



(4.8) 



As expected, all amplitudes vanish in the soft-pion limit. Moreover, resonance contributions 
vanish at zero recoil, as predicted by the heavy mass limit. 

C. 5 — > D-kIu decay rate 



In the analysis of i?^4 decays it is convenient to use the momentum combinations 



P =PD+Pn, 

Q^Pd -Ptt, 

L=Pe+ Pu: 

N ^pe-pu- 



Pd = Mdv', p^ = p, 



(4.9) 



In terms of these variables, the most general form of fl^ is 



19 



-He.^.^L'-Qpp'' + FP. + GQ. + RL, 



(4.10) 



where H, F, G and R are form factors dependent on the three invariants u, p-v and p-v'. 
These form factors are easily obtained from the exphcit expressions of the non-resonant and 



resonant parts of the amphtude given in eqns. ( [4 .21 ) and ( [4.4D . The exphcit expressions 
for the form factors are given in Appendix C. The assumption that the leptonic current is 
conserved imphes that the term proportional to R does not contribute and can be ignored. 

The squared modulus of the decay amplitude, after summing over the lepton polarizations 
and neglecting higher order terms in the pion mass squared is given by 

Y: \T\' = k' {\F\'[\{MISd.,S,)-4{N-PY_ 

+ IGp [4 {L-QY + 4 SeiSn. - 2M|,) - 4 {N-Qf 



spins 



-K St - - 5,(2M^ - 5z5.)(M^ ~S,- Sd. 



+ 1^1' 

-{L-Qf S,Sd. + L-QMlS,{Ml-S,- Sd.) 
+ SjSoA'^Ml - Sd.) - {e.^p.L^^N-ppQ'^)'' 
+ 4 Re{FG*) [-2 Ml S^ + L-Q{Ml - Se - Sd.) - 2 N-PN-Q 
+ Re{FH*) [4MlSiN-P -2L-Q N-P{Ml- Si - Sd.) 
+ N-QX{Ml,Sn.,Se) 

+ Re{GH*) [4SiN-P{2Ml- Sd.) -^{L-Q)^N-P -4MlSeN-Q 
+ 2L-QN-Q{MI-S,-Sd.) 
+ A\m{2F*G + F*HN-P + G*HN-Q) t^^p^L^'N^pPQ"} (4.11) 

The invariants appearing in this expression are 

P^ = Sd-k, 

= -N^ = Sr, 



P-Q = Ml -Ml 

P.L = ^{MI-S,-Sd.), 
L-N = 0. 



(4.12) 
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and A is Kallen's function. In order to obtain explicit expressions for the remaining invari- 
ants, namely , P-N, N-Q and L-Q, it is convenient to use as independent variables 
the quantities Sd-k, Se, and the angles 9.,,, 9e and (p. From eqn. ( [4.12| ), Sd-k and Se are the 
invariant mass squared of the nD and iu pairs, respectively. The angles 6,^, 6e and are 
illustrated in figure |^. 6'^ is the angle between the pion momentum and the direction of ? 
in the cm. frame of the ttD pair, $£ is the angle between the lepton momentum and the 
direction of L in the cm. frame of the iu pair, and is the angle between the two decay 
planes defined by the pairs (j}^,Pd) and (j)e,pt) in the rest frame of the 5-meson. This is 
the set of variables initially introduced by Cabibbo and Maksymowicz in the analysis 
of Ki4 decays. 
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FIGURES 




FIG. 1. Kinematic variables and angles. 
The remaining invariants, are 

+ cose,(M| + Sd^ - Se) X^'\Ml, Sn^, Ml) 

+ COS X^/^Ml, Sd., Se) \^'\mI, Sd., Ml) 
^ sm^e^\{Ml,SD^,Ml), 



P-7V = ^ cos^, \^/\Ml, Sd., Se), 

L-Q = ^ [mUm'b - Sd. - Se) 

+ cos^, \^'\Ml,Sn.,Ml) X^'\Ml , Sdtti Se) 
^■Q=^ cos^, X'/\Ml, Sd., Se) 

+ COS ee cos 9^ X'/\Ml, Sn., M^ 

X (SeSo. + - Sj - Si, + A(M|, 5^,,, 5,)) 
cos (/) sin Oe sin A^/2(M|„ Sd^, Ml), 



Se 



S 



Dtt 



t,,,.P''Q''LPN- = -o t/TT^ >^^'\Ml, Sd., Ml) X'/\Ml Sd., Se) 

Z V ODtt 

X sin^sin^^sin^jr- (4-13) 
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Since the form factors depend on only one of the angles, namely ^^r, in the expression for the 
partial width of the decay the integrations over the angles (j) and On can be performed 
explicitly. Following standard steps, the differential partial width of interest can be expressed 
as 



d cos d-^dSD-wdS I 2Mb 



J^Sct^jSi) / d(f) dcos6i 
Jo J -I 



where 



X \TWSD.,S,,e^,ee,(t)), (4.14) 

spins 



2, for charged pions 

K = <( , (4.15) 
1, for neutral pions 

and the Jacobian J{x, y) is 

J{x,y) = ^uj^y^2 A'/'(M|,x,y) \^/\Ml,Mlx) \^/\Q,Q,y). (4.16) 

For our purposes the interesting differential partial rate is dV b^^I dSD-n which results from 
integrating eqn. ( |4.14]) over 5*^ and 6'^ with no kinematic cuts. 



Y). B ^ D*Titu decay rate 

The tensor fi^^ can be expressed in terms of twelve form factors as 
% % % 

+ F,P^{P - Q), + F^Q^iP - Q), + F^P^L, + F^Q^L, + Kg^ 

^ A. 5 ^ A. 5 

+ ^p-SpcrP Q^L'^{P — Q)u + ^pSpaP Q^L'^Li, 

+ '-Gf e.spaP'Q'L'^Pp + ^Gf e^SpaP'Q'L^Q,. (4.17) 

In writing this form, we have neglected terms that vanish upon contraction with the con- 
served leptonic current and with the D* polarization vector e^, . The explicit expressions 
for the form factors resulting from eqn. ( |4.17| ) are given in Appendix C. 

It is straightforward to calculate the modulus squared of the resulting decay amplitude 
summed over the polarizations of the D*. Since the result is lengthy we prefer not to display 
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it here. The partial width is given by an expression similar to that of eqn. ( [4.14| ) with the 
appropriate replacement of the squared amplitude. 



V. RESULTS, DISCUSSION AND CONCLUSIONS 

If we throw caution to the wind and apply our calculation to all of phase space, we find 
that the decay rate for B DttIu ranges from 1.09 x 10"^^ GeV to 1.13 x 10"^^ GeV. The 
upper limit corresponds to including all the multiplets in the calculation, while the lower 
limit arises from including only the (0~,1~) and (0"^, 1"^) multiplets. These decay rates 
correspond to branching fractions of 2.1%, somewhat larger than, but largely in agreement 
with the analysis of Cheng and collaborators ^ . We see, therefore, that the inclusion of the 
higher multiplets does not profoundly affect the total decay rate of S — Dniu. The effects 
on the spectrum, and on the decay B D*7Tii' are somewhat more striking, however. 

Performing the same integration for the decay B D*7r£z/, we find that the total rate 
varies between 2.7 x IQ-^^ GeV {BR = 5.0 x IQ-^) and 1.7 x lO^^^ GeV {BR = 0.3%). 
Thus, inclusion of the resonances makes a significant change to this decay rate, increasing 
it by a factor of about 6. In both cases {B — > DTciu and B D*Tilv), the change in the 
quark model parameters from set I to set II has little effect on the integrated decay rates, 
but makes significant differences to the spectra obtained. 

In a series of figures we show the differential decay rates into a charged pion q^^* 
as a function of Sd-k-, for values of Sd-k between the threshold of {M^^^ +m7r)^ and 10 GeV^. 
This covers a bit beyond the domain where the soft pion limit may be safely applied; for 

~ 0.5 GeV, the soft pion limit holds up to Sdt, ~ 6.5 GeV^. In each of these figures, we 
show the spectra resulting from both sets of parameters, with set I corresponding to figures 
N(a), and set II to figures N(b). In addition, we normalize by dividing the differential decay 
width by the total semileptonic decay width B — ^ Dlv, calculated in the same model. In 
this way, we can lessen the impact of model dependences that enter through form factors 
and coupling constants. We note, however, that using \Vcb\ = 0.043, we obtain branching 
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fractions of 1.8% (1.7%) (here, and in all that follows, the first number is obtained using 
the parameters of set I, while the second number, in parantheses, is obtained using the 
parameters of set II) for B —>■ Diu and 4.6% (4.6%) for B —>■ D*lv, in surprisingly good 
agreement with experiment. We have not compared our results with the differential decay 
rates for these decays, however. 
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FIG. 2. Y{B^D^v) dSo^ ^ function of Sdtt for B Dir^eu. The different curves correspond 
to different combinations of resonances, as explained in the figure. 





FIG. 3. Y{B^Dev) dSo^ ^ function of Sd-k for B D*TT^eiy. 
Figures ^ and |^ show the spectra for the decays B Dniu and B D*nii', respec- 
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tively. Each of these figures shows four curves, corresponding to the inclusion of various 
combinations of multiplets in the analysis. One very interesting feature in these spectra is 
the narrow structure due to the (0~, D meson multiplet. We will discuss this feature in 
some more detail later. Another notable feature is the depletion in the differential and total 
widths ol B ^ Dniu that arises from the interference between the (0~, 1~) and (0"*", l"*") 
multiplets. In the case of the decay B — > D*TTii', the depletion caused by this interference 
at values of Sd*-k < 8.5 GeV^ is more than compensated by the enhancement that occurs 
throughout the rest of phase space. This latter enhancement should however be taken with 
caution, as it lies beyond the range of applicability of our approximation. 




FIG. 4. The effect on the decay B Dir^eu, of changing the couplings and total widths of 
the states in the (0+, 1+) multiplet. 
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FIG. 5. The effect on tlie decay B D*tt eu, of changing the couplings and total widths of 
the states in the (0+, 1+) multiplet. 

In figures |^ to ^ we examine the effects of the values of some of the parameters on 
the spectra. Figures |^ and ^ show the effect of changing the coupling constant and width 
of the states in the (0"^, l'^) multiplet. The values we have obtained in our model are 
ai = —1.43 (—1.22) and T = 1.04 (0.756) GeV. This width may seem very large, but it 
is at least consistent with other model calculations. Nevertheless, we have investigated the 
effect of using smaller widths, namely 500 MeV and 200 MeV. ai is changed to -0.99 and 
-0.63, respectively, in keeping with these changes. We see that the effect on the spectrum of 
B — > Diriu is quite striking, especially at the narrower of the two widths. The effect on the 
spectrum of B ^ D*Tilv is similarly striking. The total width of the B D*TiC.v decay is 
strongly affected by these changes, changing from 1.7 (1.7) x 10^^^ GeV to 1.0 (1.2) x 10^^^ 
GeV and 4.9 (5.9) x 10"^^ GeV as the total width of this multiplet changes from 1.04 (0.756) 
GeV to 500 MeV to 200 MeV. In comparison, the total width of i? ^ Dniu is essentially 
unaffected by these changes. 
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FIG. 6. The effect on tlie decay B D-K^eu, of changing the couphngs and total widths of 
the states in the (0^, 1^)' multiplet. 





FIG. 7. The effect on the decay B D*Tr^ei', of changing the couplings and total widths of 
the states in the (0", 1^)' multiplet. 



In figures |^ and ^ we illustrate the effect of changing the total width of the (0~, 
multiplet from 191 (174) MeV to 130 MeV, accompanied by a change in from 0.69 (0.66) 
to 0.57. We see that the narrower state would provide a clearer signal for experimentalists. 
The possibility of observing this pair of states in Bi^ decays will depend strongly on the 
value of the total width and on as, as well as on the effects that various experimental cuts 
will have on the spectra we illustrate. While we do not study the effects of cuts in this work, 
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we can estimate the number of events that one may see in the proposed S-factory. 

The integrated width under the peak (from about 5.33 to 6.0 GeV^) is 2.7 (3.9) x 10^^^ 
GeV, corresponding to a branching fraction of 5.2 (7.4) x 10~^, or about 5.2 (7.4) x 10^ 
events, assuming production of 10® B mesons. Subtracting the width that corresponds to 
a 'smooth' background leaves 220 (7000) events in the peak alone. A similar exercise in 
the case of the D*7r£i/ spectrum yields a width of 7.4 (10.3) x 10^^'' GeV, corresponding to 
14200 (19900) events. Removing the smooth background leaves 1800 (4300) events in the 
resonant peak. 

While the dependence on the model parameters is clear, these estimates suggest that a 
study of the spectra of B ^ DTriu and B D*Tilv offer some opportunity for discovery 
or confirmation of the resonances of the (0^, 1^)' multiplet. Note that if we include the 
expected small mass difference between these two states, the single peak in these figures will 
become two peaks that are very close together (separated by about 0.32 GeV^ if the mass 
splitting is about 60 MeV). The net effect would be a broadening of the structure that we 
have in our spectra. 

In conclusion, we have studied B^^^ decays in the soft-pion limit using chiral perturbation 
theory and heavy quark symmetry. The resonances which give leading contributions in 
this limit have been included, and shown to be important in determining both the rate 
and the shape of the spectrum. The narrow (0~, 1~)' resonances show up as a peak in the 
Sot^ spectrum, and will likely be difficult to observe in B ^ Diriu. The possibility for 
detection or confirmation in B D*t[Iv is more promising. The wider resonances of the 
(0"^, 1"^) multiplet show some effect on the total rate, but are not likely to be identified 
from the spectrum. Preliminary indications are that they may be identified at Aleph using 
the topology of the B^^ decays |2T|. The effect of these broad resonances is much more 
pronounced for B —>■ D*Ttlv than for B Dniu, although the effects on both spectra are 
quite clear. 
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APPENDIX A: STRONG INTERACTION TRANSITION VERTICES WITH A 

SINGLE PION 



1- i p» 



1- i P» 



3 



FIG. 8. Vertices obtained from the chiral Lagrangian for soft neutral pion emission. Vertices 
for charged pion emission are a factor \/2 of those shown in this figure. 

In this appendix, we give the explicit expressions for the vertices where one pion emission 
takes place according to Similar results hold for B mesons. The vertices are shown in 
figure |. 
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APPENDIX B: THE CHARGED CURRENTS 



In this appendix the explicit expressions for the charged currents displayed in eq. (p.l2|) 
are presented. They are 
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1^ 



D\v')iv + v'),B{v), 



X 



X 



0- ') 
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gup{v + v')^ - g^pV^ - gf^yv'p + it^apy{v + v'Y B*P{v), 
-p,{v) Dl{v'){v-v'),B{v), 

pi{u) D*^'"{v') [ie^^a/BV^'v''^ + 5(^j.(z/ - 1) - v'^v^] B{v), 
pi{v) Diiv') [-te^a/SpV'v'^ - g^piu - 1) + v^v'^] B*p{v), 

X 'tep.apVv'f'iu - 1) + gp,{i^^ - 1) - v,{{2 + v)v'^ - ^v,)] B{v), 

X -v^g^p{iy-l) + VuVpv'p + iep^apVpv"v''^ B*''{v), 
= ^p,{u) Dl-'^v') 

X [(^ + v')pgya{v - 1) - "iVyVpv'^ + 2v^gp^{u - 1) - gpyv'„{v - 1) 
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■e\v)D'\v'){v + v')pB{v), 
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B*P(v) 



(Bl) 



The effective currents where a S-meson resonance decays into a ground state D-meson 
are easily obtained simply taking the hermitian (minus the hermitian) conjugate of the 
vector (axial-vector) portion of the currents displayed above followed by the interchange of 
symbols B ^ D and v ^ v' . 

APPENDIX C: THE FORM FACTORS 



In this appendix we give the form factors needed in eqns. ( [4.10D and (|4.17| ). The 
non-resonant and the resonant contributions are displayed separately. 



1. B ^ DttIv 



If we write 



= -ih MbMd e^.p^v'^v'Ppl + Ai p^p + A2 Mbv^ + A3 Mdv' 



(CI) 



the form factors in equation (^.101) are 



H = -K 



F = A2 + -(Ai+A3) 
^ = ^(^3-^0, 

R = Ao. 



(C2) 



From eqn. ( [4 .21 ) we obtain the non-resonant contributions to the form factors as 

9 ^(^) /I 1 \ 



h 



Ai 



NR 



NR 



2Fq MbMb, \Pt,-v + SrriB — ie Pt,-v' — Smj:, + ie 



9 
2Fo 



(1 + ^) 



.p-K-v + 5mB — Pt,-v' — 6mD + i€ 
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. 9 ^{J^) ( Pn-V+Pn-V' 

^2NR 



AsNR - -TTw TT- z—:, — Tz — 7— ■ y^'^) 



2Fq Mb \Pn-v + SniB — ie ^ 
9 ^{i^) ( Pn-v + Ptt-v' 

These results are the same as those obtained by Lee and collaborators J^, and by Cheng 
and collaborators [|8|. 




FIG. 9. Non-resonant (a) and resonant (b) diagrams contributing to i? ^ DttIu. The dashed 
line represents the pion and the dotted line which emerges from the electroweak charged current 
carries the momentum of the pair. Horizontal solid lines correspond to four velocity v and the 
oblique ones to v' . 
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From eqn. 



the resonant contributions are 
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(C4) 



2. B ^ D*tt£u 

The most general form for the tensor in terms of the vectors v^, v'^ and p.,^^ is (terms 
which vanish upon contraction with t*Q{y') are not displayed) 

n^,{v,v\p^) = t^,p, [/ii MbMovPv"' + h2 MBvPpl + /i3 M,5t;'X] 

+ /i MbV^^P^u + /2 Mdv'^P^u + /3 PttmPtti. + /4 M^V^V^ 

+ /g MbMdv'^v^ + fe MbP^^^v^ + k g^^ 
+ ^ ef,spa v^v'^pI {gi p^u + g2 MbVu) 
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% / \ 

+ ^ ^v&pa V^v'Ppl (5(3 MbV^ + 5-4 Mdv'^^ + 5(5 p^^J . 



The form factors appearing in eqn. ( [4.17| ) are related to the ones in this expression 
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"2 {hi + h2), 
I i-hi + h), 

1 (/l + /2 + ^/3 + /4 + ^/5 + ^/e 

^ (/2 — /s + /s — /e) , 

2 if 5 ~ fe), 



1 



AMbMd 
1 



(^1 + ^2), 



-^2, 



(^3 + ^(^4 + ^5)), 
(^4-^5)- 



1 



AMbMd 
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(b) 



FIG. 10. Non-resonant (a) and resonant (b) diagrams contributing to B 



From eqn. ( |4.6|) the non-resonant contributions to the form factors are 



hlNR 
h2NR 
h^NR 
flNR 



FqMbMd p^-v + SniB - ie'' 



FqMb Pt,-v + 5mB - ie 



FqMd \Pn-v + SmB-ie p^^-v' + it^ 
9^^) ( 1 1 



2FqMb xPtt-v + 6mB — ie Pt,-v' + drrin + 



hNR — -TT~ /lA^-R' 
fsNR = fiNR = 0, 
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^'"^"^ 2FoMbMd V Pn-v + 6mB-te) 

f 9^^) ( 1 1 

Je NR 



2FqMb \Pt,-v + 5mB — it p^-v' + ie J 

2Fo \Pt^-v + Sttib — ie p-^-v' + ie 

91 NR = 0, 

92 NR = 0, 
92, NR = 0, 



94 NR 



FqMd Pn-v' + ie 

96NR = 0. (C7) 

These results are the same as those obtained by other authors 

Finally, the resonant contributions are obtained from eqn. ( and are 



h 



IR. 



h2R 

hsR 



FqMbMd \—Pt,-v — 6mi pT,-v' — 6miJ 
SFqMbMd V Pn-v + 6fh2 2 Pt,-v' - 5rfi2 

FoMBMDPn-V + 5m3' 

a2 {1 + P2{i^) 1 (y3^^^K^) 



lii; 



SFqMb -Pn-v - Sm2 FqMb {Pn-v + Srhs) ' 

3FqMd \p^-v + Sm2 2{p^-v' - 5fh2) ) 

as ^^'K'^) ( 1 \±J^\ 

FqMd \p^-v + 6m3 p^-v'-6msj' 



GFqMb \p.„-v + 5m2 Pt,-v' - 5m2^ 

as e'K^) ( 1 1 ^ 

2FqMb \Pt,-v + 5m-i p^-v' - 6m 



3, 



_Mb 
J2R — -ry- JiR, 
Mb 

/siJ = /4_R = 0, 

ai pi(z/) / p^-f p^-f' 



2FoMbMd \Pn-v + 5mi p-„-v' — Srhi 
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2FqMbMd V Pn-v + 5m2 p^^-v' - 5m2 



f&R — 
+ 



2FoMbMd ypTT-v + Srhs p^^-v' - Sm^ J ' 



as e'Kiy) ( 1 1 \ 

2FqMb \p-„-v + 5rh^ p^r-v' — Srhs J ' 

Q!i pi(z/) (z/ - 1) / p^-v p^-v' 
f^R^ 777^ + 



2-fo \p-„-v + 6mi pT^-v' — 6miJ 



_ Q2 P2('^) 1 

^'"^^ 2FoMbP^-v'- 67712' 

93 R — —g2R, 



+ 



SFqMd yP7r'f + <^"^2 2 Pt^-v' — 61712^ 



5'5ii = 0. 
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